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Abstract. - We study the T — normal to superconducting transition of Josephson arrays 
with on-site dissipation. A perturbative renormalization group solution is given. Like the 
previously studied case of bond dissipation (BD), this is a "floating" to coupled (FC) phase 
transition. Unlike the BD transition, at which only nearest-neighbor couplings become relevant, 
here all inter-grain couplings, out to infinitely large distances, do so simultaneously. We predict, 
for the first time in an FC transition, a diverging spatial correlation length. Our results show 
the robustness of floating phases in dissipative quantum systems. 



Coupling dissipation to quantum systems has interesting consequences. Quite generally, 
dissipation suppresses quantum fluctuations. This is evident at the level of a single macro- 
scopic quantum 'particle' whose tunnelling probability out of a metastable state is suppressed 
by dissipation [1]. For the quantum particle in a double well [2,3], or in a periodic [4,5] 
potential, Ohmic (i.e., linear) dissipation can suppress quantum fluctuations enough to cause 
a dissipative 'quantum to classical' phase transition. In the classical phase, all quantum fluc- 
tuations are quenched and the system is spontaneously trapped into only one of the potential 
minima. 

From the point of view of understanding effects of dissipation on extended systems, a 
natural question to ask is what happens when such zero-dimensional dissipative systems are 
spatially coupled in a lattice. Understanding such effects has important practical applications. 
Such dissipative effects are thought to be at the heart of the physics of granular supercon- 
ductivity [6]. They are also important in the context of decoherence in a qubit, which, it 
has been proposed [7], can be realized in a Josephson junction array . If the coupling of 
dissipation remains local, then it is natural to expect that quantum fluctuations can also be 
locally quenched by dissipation. At zero temperature, coupling such classical systems spatially 
should give rise to an ordered state of the extended system. This dissipative quantum phase 
transition (QPT), however, is expected to retain some of its local character. 
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Early work on these questions [6,8-13] was recently extended in Ref. [14] in the context of 
a Josephson array with dissipation coupled to the bond phases, that is, the order parameter 
phase differences of two consecutive superconducting grains; see also Ref. [15]. This situation 
can be realized when the junctions are shunted by Ohmic resistors. It was shown that the metal 
to superconductor dissipation-tuned QPT occurs via a floating to coupled phase transition. 
Interesting local behavior in the bond dissipation model has also been addressed in Ref. [16]. 

The case where dissipation couples to the site variables - that is, the phases of the su- 
perconducting islands - directly has also been addressed in the literature [17-19]. Such local 
damping of the phase happens [17] when Cooper pairs can decay into localized pools of normal 
electrons in the substrate. It can also be built into an artificial array with each island shunted 
separately to the ground [19]. Such a 'shunt to the ground' also occurs naturally in many 
granular systems. The model has been mapped on to effective LZW field theories [17, 19] to 
study the normal to superconductor phase transition. 

Working directly with the order parameter phase variables, we give in this paper a sec- 
ond order perturbative renormalization group (RG) solution for the site dissipation problem. 
We show that a floating normal metallic phase, which is even more striking than in the 
bond-dissipation case, and a floating to coupled QPT occur in this case. In contrast to 
bond-dissipation, here all couplings between grains become relevant simultaneously and the 
transition is controlled by a fixed, infinite-dimensional critical hypersurface with continuously 
varying critical exponents whose dependence on dissipation we calculate near the threshold 
dissipation. We also predict, for the first time in a floating to coupled phase transition (three 
other examples of which are known [14,20,21]), a diverging spatial correlation length. The 
RG flows and phase diagram we derive are displayed in Fig. ^ in the V — a plane, where V is 
the nearest neighbor Josephson coupling and a is the coefficient of dissipation. 




Fig. 1 - Phase diagram and perturbative RG flows in the V — a plane in any dimension D > 0. 
N and SC indicate normal and superconducting phases, respectively. The normal metallic phase is 
the floating phase. The thick dashed lines correspond to lines of fixed points. The cross indicates a 
special point on the a = axis denoting a (D + l)-dimensional XY transition. The part of the a — 
axis from the cross to infinity is a line of fixed points in D — 1. 



Our results demonstrate the model-independence of floating phases in dissipative Joseph- 
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son arrays. They also indicate that dissipative quantum systems may be the best places to 
experimentally look for such phases, which may otherwise be very difficult to obtain [22]. 
Moreover, they provide a rare example in statistical mechanics of a critical fixed locus that 
exists in any dimension; there is no upper critical dimension in the problem. 

On-site dissipation coupled to XY spins has also been treated in Refs. [23,24]. There, 
dissipation is coupled locally to a </> 4 -theory, which retains the compactness between and 
2tt of the phases of the order parameter. The problem is then studied using an e-expansion 
around the upper critical dimension. In our model, the compactness is explicitly broken, as is 
expected in Ohmic dissipation [25]. We believe the difference of our results from Refs. [23,24] 
is due to this different treatment of the phase variable. 

The action for the Josephson junction array coupled to on site Ohmic dissipation reads 

S / h = EE^ + ^n\)\0(r^n)\ 2 
r n 

+ V /dr[l-cosA%^(r)]. (f) 

(r,r f ) 



Here C is the capacitance of the grains arranged on an arbitrary .D-dimensional lattice 
with lattice points r. 9(r,r) is the phase of the superconducting order parameter on the 
grain at r at imaginary time r, and 9{f,uj n ) is its temporal Fourier transform. The bosonic 
Matsubara frequency u) n — 2im/ [3 with n an integer, and /3 the inverse-temperature, (r, r*) 
denotes nearest neighbor pairs, and A9^p = 9{r) — 9(r'). The term containing a arises 
from integrating out a harmonic-oscillator bath with the linear spectral-function required to 
produce Ohmic resistance [f]. In r-space, this term reads £ J dr J dr' ^^P^- f, which 
is a long ranged interaction in time. We restrict the sum on the Matsubara frequencies to 
those satisfying \uj n \ < lu c , where lo c is an ultraviolet cutoff . 

To study the dissipation-tuned quantum phase transition in this model, we use a renormal- 
ization group that is perturbative in the spatial coupling V [6, 14]. This RG procedure starts 
by dividing the field 9(r,r) into slow and fast components # s (r, r) and 9f(r,r), respectively, 
defined by having the temporal Fourier transform of 9f{r 1 r) equal that of 9(r, t) for "high" 
Matsubara frequencies obeying ^ < \uj n \ < oj Ci and equaling zero for all smaller Matsubara 
frequencies, while the temporal Fourier transform of 9 s (r, r) equals that of 9(f,r) for "low" 
Matsubara frequencies obeying |w„| < We can then write the partition function Z as 



Y[ve s ( 



r,r) exp 



s^/n + He- s ' /n )o f ]. (2) 



Here Zq is a normalization constant, Sq is the slow-frequency component of the quadratic 
part of the action, S' contains the spatial coupling term V , and (. . .)o/ denotes averages over 
the fast components. After computing the averages, we rescale r, r' = r/b, b being the scale 
factor, to restore the original frequency cut-off. Finally we redefine the coupling constants to 
complete the renormalization. Note that we do not rescale the fields 9(r, r); this is to avoid 
introducing a factor multiplying A9?,p (r) in the argument of the cosine in Eq. Q . With this 
choice to not rescale 9(f, r), the dissipation term a in Eq. is dimensionless. Furthermore, 
because it is non-analytic in u> n , it gets no "graphical" corrections (i.e., it is unchanged by 
the step of integrating out the "fast" fields). Hence, it is held fixed by the RG. The term 
associated with C has r-dimension —1; it is irrelevant in the RG sense. The only place where 
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C appears is in the frequency integrals such as, 

duj 1 = 1 a(3 

Thus, ^ acts as an upper frequency cut off, and we associate it with uj c , ui c = 
Writing b = e l , where I > is infinitesimal, we have, to one loop order [6, 14]: 

f-d-i)V (4) 

This rcursion relation implies that the a-axis (V = 0) is a line of stable fixed points for 
a < 1 and a line of unstable fixed points for a > 1. Defining a c = 1, for a < a c , the 
barrier V between the different potential minima of the cosine potential is irrelevant, and 
so the fields (9(r, r) on each site fluctuate quantum mechanically independently of those on 
all the other sites. For a > a c , on the other hand, the barrier V grows, quenching these 
fluctuations, and the system becomes a global superconductor. Therefore, at a = 1, there is 
a metal-to-superconductor quantum phase transition. 

Before we go to the higher order corrections to Eq. 0] and discuss the phase diagram, we 
first show that on the metallic side of a = 1, the local phases 9(r, r), which are critical in the 
imaginary time variable, exist in a so called 'floating' phase. The floating phase for the array 
is defined as a phase where all spatial Josephson couplings, including the longer ranged ones 
that nominally couple spatially far separated grains, are irrelevant [14,20]. From Eq. ^with 
irrelevant V, the average of the order parameter on a grain is zero at T = 0, (exp(i6>(r, r))) = 
1/(^/3)2^. The on site unequal time correlation function decays algebraically, 

(expi(0(f, r) - 9(r, 0))) = — ^ (5) 

(UJ C T) a 

Now, a generic arbitrary range Josephson coupling among pairs of grains takes the form: 

Sj = -fdTj2 J{?) cos(6(r, T )-e(r + r*,T)), (6) 

r,r' 

where \r'\ gives the spatial range of the coupling. When V is irrelevant, the average of Sj 
scales as (Sj) = /3(exp(i0(r, r))) 2 - P 1 "'. At zero temperature, this term is irrelevant for 
a < 1. Thus, all spatial couplings, including the nearest neighbor coupling V, are irrelevant 
throughout the metallic phase. Interestingly, they all become relevant simultaneously at a = 1. 
Hence, for a < 1, the system is in a floating critical phase. 

The first order flow equation, Eq.0J only establishes a phase transition. To know the phase 
diagram and the critical properties, such as correlation length and time exponents, we have to 
go beyond first order. To obtain higher order corrections, we perform a cumulant expansion 
of Eq. in second order in V. By the method outlined in Ref. [14], we find higher order 
contributions to Eq. 01 coming from the closed loops on the lattice. At higher orders in RG, 
longer ranged pairwise interactions J(r) of Eq. are successively generated. In a consistent 
treatment of the RG, all these terms should be included in the starting Hamiltonian. Since 
the quadratic part of Eq. ^ decouples in space, it is easy to see that all the J(r)'s have the 
same linear order recursion relation as that of the nearest neighbor interaction V, Eq. At 
second order cumulant expansion of (J2J, J(r) for a particular r is generated by all possible 
combinations of J(r') and a J(f — f") for which r, f* and r — r 1 lie on a triangle on the 
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lattice [14]. Treating V as a special case of J(r) (that corresponding to the nearest neighbor 
vectors d\, Jifli) — V), we get, up to quadratic order in the couplings, 

' iJlF] + (7) 



dl 



Here A is a positive number, A ~ f Q d{r — r')[exp(— G? (t — r')) — 1], where G s (j — r') =< 
Of(Q, r)9f(0, t 1 ) >of is the unequal time correlation function of the fast modes [14]. To avoid 
generation of spurious long ranged interactions, a smooth cut off prescription is adopted to 
calculate G^(t — t'), which then makes it exponentially decaying with |r — t'\ [5,26]. Since A 
depends on this cut off scheme, physical quantities such as exponents should be independent 
of A. 

It is essential for our subsequent analysis that all of the factors J(f)J{f— f) enter the 
recursion relation for J(r) with the same coefficient A. This is a simple consequence of the fact 
that, in contrast to the bond dissipation case, here, in the absence of the Josephson coupling 
terms, fluctuations of 0{r, r) on any site r are completely independent of those on all other 
sites f . As a result, the averages over 9(f, r)'s on different sites that must be computed in 
the perturbation theory to obtain A are the same, regardless of which sites r and r — r 1 are 
being considered. 

Because of the fact that A is a constant that factors out of the sum on f* in Q, that 
formidable looking set of infinitely many coupled non-linear recursion relations decouple in 
Fourier space. Defining the Fourier transform J(q) of J(f) via J(r) = J BZ d d qj(q) exp(i<f ■ r) 
where the J BZ d d q is over the Brillouin zone of the lattice, and defining e = ^ — 1, Eq. [7| 
reduces to 

d J^L = -eJ{d) + AJ*{d). (8) 

Let's consider the metallic side of a c = 1, where e is positive. Equation[S]has a stable fixed 
point at J(q) = 0, and a critical fixed point at J(q) = e/A. Fig. [21shows the flow of J(q). 



m 



— >>><<< — « x » 

e/ A 

Fig. 2 - RG flows for J(q) as implied by Eq. El 



It follows that as long as the maximum of J(q) in g-space is less than V c (e) — e/A, 
the system is in the floating metallic phase. As shown in Fig. |3| when the maximum lies 
above V c (e), J(q)'s for a band of q values become relevant and the system becomes a global 
superconductor. To illustrate this, let's consider a particular form of the bare J(r), namely, 
that for nearest neighbor interactions only, on a D-dimcnsional hypercubic lattice, J(r, I = 
0) = VJ2a 5(r ~ (Ti), where the d^'s are the set of nearest neighbor vectors on this lattice. 
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This gives a bare J(q) of the form, 

J(q, l = 0) = Vj2cos(q- CLi). (9) 

di 




Fig. 3 - A typical form of the bare J(q) in g-space. The dotted horizontal line indicates V c (e). When 
the maximum of J(q) at q — lies below the dotted line (solid curve), the system is in the floating 
metallic phase. When the maximum lies above the dotted line (dashed curve), J(q)'s in a band of q 
values, \q\ < q c , become relevant, and the system becomes a global superconductor. 

As shown in Fig. [31 as long as zV is less than V c (e), where z = 2d is the coordination 
number of the lattice, all J(q)'s are irrelevant. When zV is greater than V c (e), J(q)'s in a 
band of q values are relevant. We expand J(q, I = 0) in the neighborhood of q = 0, 

J(£/ = 0)~ J(q = 0)-B\q\ 2 , (10) 

where J(q = 0) = zV and B — 2Va 2 (with a the length of a nearest neighbor vector) 
in the current example. For more general lattice types, and further than nearest neighbor 
interactions, the form of Eq.^Jwill still hold for small q, provided only that the bare Josephson 
couplings J(r) are such as to produce a net ferromagnetic interaction (by which we mean one 
that favors a spatially uniform #(r, r)), and short-ranged. Indeed, that J(q) has its maximum 
at q = , and that that maximum is a simple quadratic maximum as in Eg. 1101 can be taken as 
a definition of ferromagnetic, short-ranged interactions. Furthermore, it will be satisfied by any 
system that has, in the bare Hamiltonian, only nearest-neighbor ferromagnetic interactions, 
whatever the lattice type. 

Equating J(q, I — 0) in the neighborhood of its maximum with V c (e), we get the edges of 
the band of q values \q\ < q c , q c — — y==(V — Vc ^ ) ^ ; for which the J(g)'s are relevant. This 

defines a diverging correlation length £ ~ q^ 1 ~ (V — Vc ^ )~* with a mean field like exponent 
\. We have thus established, we believe for the first time, a diverging spatial correlation length 
at the floating to coupled (classical or quantum) phase transition. This length £ could be 
measured experimentally using the AC Josephson effect. Specifically, we predict that for a ID 
chain of N junctions, with opposite ends held at a small temporally constant voltage difference 
AV , an oscillating Josephson current Ij(t) will occur, whose frequency Q is independent of 
N for << £, while being inversely proportional to N for A^ >> £. The crossover between 
these two behaviors determines £. The linear form of V c (e) = 4, where e = a c — a, determines 
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the RG flows and the perturbative phase diagram in the V — a plane as shown in Fig. ^ 
Finally, linearizing Eq.[Hj around the critical fixed point, we get the correlation time exponent 
v T , £ T ~ \V — 14(e)| _I/T , with the continuously varying critical exponent v T — - = - 1 _ a . This 
defines a dynamic exponent z, defined by £ r ~ £ z , to be z = 2v T = - 2 _ a . 

To reiterate the differences of the present calculation from the case of bond dissipation [14], 
we note that although the linear order RG equation, Eq. 0] is identical in both models, the 
second order equations, Eqs. Q|SJ, are different. Because of the constancy of the number A, 
we are able here to solve the set of infinitely many coupled non-linear equations exactly by 
going to the Fourier space. This allows us to compute important critical exponents such as v 
and v T . We are also able here to demonstrate the divergence of the spatial correlation length, 
to our knowledge for the first time, in an FC transition. Note that, in the corresponding 
problem of the coupled XY planes [20], or the case of bond dissipation [14], this was an 
interesting open problem because of technical difficulties. 

In systems where the 'shunt to the ground' is realized, either by a ground plane, or in- 
trinsically in granular systems [17, 19], the FC transition is tunable by varying the shunt 
resistance (which varies a in our model). By measuring the value of a at which the transition 
happens for a variety of different V's, it should be possible to test our prediction for a c , and 
that the phase boundary comes into the a-axis at a c with a finite, non-zero slope. A similar 
phase diagram has been obtained experimentally in a dissipative array in Ref. [27], but in 
that experiment a bond dissipation model was more applicable. In systems where the present 
model is applicable, the non-universal exponent v T can be directly measured by measuring 
a characteristic temperature, T ch cx \V — V c {e)\ VT , in either phase close to the phase bound- 
ary. In the superconducting side T c h can be the superconducting T c , while in the resistive 
side, it can be the temperature where the resistance shows a minimum [8]. On the resistive, 
metallic side of the transition, the floating character can be established by measuring the 
current-voltage characteristics which are expected to be a non-uniform power law depending 
on a [28]. Finally, although the correct model for capacitances of a real experimental system 
may also include longer ranged capacitors, in contrast to the on-site capacitance C used here, 
we can show that this does not change any of our universal results. 

We thank Sudip Chakravarty for useful discussions; JT also thanks the Aspen Center for 
Physics for their hospitality while a portion of this work was being completed. This work 
was supported by the NSF under grant Nos. DMR-01-32555 and DMR-01-32726, and by 
ARO-ARDA under grant No. W911NF-04-1-0236. 
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